Simple model for reverse buoyancy in a vibrated granular system 
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Large objects, immersed in a homogeneous granular system, migrate when subjected to vibrations. 
Under certain conditions large heavy objects rise and similar light ones sink to the bottom. This 
is called reverse buoyancy. We report an experimental study of this singular behavior, for a large 
sphere immersed in a deep granular bed. A simple mechanism is proposed to describe the motion 
of a sphere, inside a vertically vibrated granular system. When reverse buoyancy is observed, the 
measured vertical velocity of the immersed object, as a function of its density, shows a simple 
behavior. With a one-dimensional mechanical model that takes into account a buoyancy force and 
the frictional drag, we obtain the rising velocity for heavy objects and the sinking rate for light ones. 
The model yields a very good qualitative and quantitative agreement with the experiment. 
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Granular matter is presently a subject of extensive 
research, due to the interesting behavior frequently en- 
countered in this kind of systems (dilatancy, segregation, 
arching, clustering, etc.), and the many applications in 
industry in areas like agriculture, food processing, geo- 
physics, pharmacology, etc., that can profit from those 
investigations [1,2]. 

Dry granular materials are difficult to mix homoge- 
neously. Industrial processing involving mixtures of dif- 
ferent grains very often lead to undesirable segregation of 
these. Properties of the grains like size, shape, density, 
etc., determine their separation in unexpected ways. It 
has become clear that segregation of grains under shak- 
ing, is a very rich but poorly understood phenomenon. 
Size segregation has caught the attention of many scien- 
tists leading to an extensive research on different aspects 
of the problem. This research has revealed the complex- 
ity involved and the many aspects of grain separation 
processes that are still unsolved. 

Shinbrot and Muzzio [3] reported that large heavy ob- 
jects, called intruders, immersed in a granular medium, 
rise when shaken, while similar objects that are light 
migrate to the bottom. This novel effect, that shows 
the complexity and richness of the behavior of granular 
matter subjected to vibrations, was called reverse buoy- 
ancy. The above authors emphasized the role of inertia 
in the upward motion of the intruders, but the cause of 
the downward motion, for light objects, has gone unex- 
plained. It has become clear that much remains to be 
learned from this intriguing behavior. 

In the migration of a large object immersed in a deep 
granular bed there may be many factors playing a role 
simultaneously. Some of the factors are: small particles 
migrating to the void left out by a large object [4,5], 



Corresponding author, e-mail address: gustav@usb.ve 



convection currents [6-9], arching [10], buoyancy force 
[11], and condensation [12]. We expect to find different 
regimes where one or more of these factors are dominant, 
depending on the way the system is shaken and on several 
characteristics of the granular medium or the intruder. 
It should be interesting to be able to separate different 
regimes and determine the relative weight of the different 
factors affecting the migration of large objects in a bed 
of small particles. This could increase significantly our 
understanding of granular matter and would certainly 
have a significant impact in industry. 

Mobius et al [13] have pointed out that the rising time 
of the intruder may be used as a sensitive probe of these 
different interactions. They also mention that increasing 
evidence indicates that differences in particle density af- 
fect size separation in granular mixtures. In this work, we 
have measured the velocity of migration for large objects 
with equal shapes, similar sizes and different densities, 
when reverse buoyancy is observed. With a very simple 
mechanical model we are able to determine the velocity 
of the intruder as a function of its density. 

For our experimental set up we used spherical glass 
beads with diameters between 210 and 250/xm, the in- 
truders were spheres of approximately 2.5cm in diam- 
eter. Hollow plastic spheres were filled with different 
materials to vary the density. The vertical vibrations 
were produced with a VTS100 electromagnetic shaker 
coupled to a Plexiglas rectangular container with vol- 
ume equal to 15 x 8 x 2.7cm 3 . The frequency used was 
/ = 11.7 'Hz, the amplitude was fixed at A = 8.5mm, 
so that the adimensional amplitude of the acceleration 
was r = A(2iTf ) 2 1 g = 4.7, where g is the acceleration of 
gravity. T was monitored using an ADXL150 accelerom- 
eter. The container was filled with the glass beads up to 
a height of 11cm and the intruder was introduced in the 
medium at a certain depth of the container. Since the 
large spheres could be seen through the Plexiglas wall, a 
digital camera was used to register the position of the in- 
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trudcr as a function of time (see figure 1). The density of 
the intruder was varied from 0.26gr/cm 3 to 7.7gr/cm . 
These correspond to a range of approximately 20% to 
500% of the density of the medium. This density was 
measured by sampling different volumes of the granu- 
lar medium and measuring the weight with an electronic 
balance. 

In figure 2 we show a typical graph of the vertical posi- 
tion of the intruder as a function of time for a heavy (a) 
and a light (b) intruder. The reference frame is located 
on the container with the positive y-axis in the upward 
direction. Figure 3 shows the velocity of the intruders 
as a function of the ratio pi/p m , where pi is the den- 
sity of the intruder and p m is the density of the granular 
medium. We can see that, for the measured densities, 
the average velocity of the heavy intruders is positive, 
and its magnitude increases monotonically with its den- 
sity. While, for the light intruders, the average velocity 
is negative and its magnitude decreases with increasing 
density. 

In an attempt to explain the above results we formu- 
late the following simple model: It is assumed that the 
medium fluidizes only in part of the cycle of oscillation. 
This was proposed by de Gennes in relation to an in- 
stability in a sand heap [14]. He proposed a series of 
passive and active regimes. A passive regime occurs in 
each oscillation, when the grains are compacted and can 
be considered as a solid, and an active regime develops 
in the same cycle, when the granular system is fluidized. 
According to this approach, during the active fraction of 
the oscillation, the system is subjected to an apparent 
gravity, reversed upward. This reversal of the apparent 
gravity increases the relative interparticle motion of the 
medium causing a fluidization. This assumption was used 
in connection with instabilities observed in shaken sand 
in a U-tube [15], a similar condition was used to describe 
a transition from a condensed state to a fluidized state 
in a column of beads [16], and to derive a criterion for 
the onset of covection [17]. When the sand is not flu- 
idized it behaves like a solid and the intruder is trapped 
and does not move relative to the granular bed. When 
it is fluidized we assume that three forces are dominant: 
the weight, the drag and the buoyancy. In the noniner- 
tial reference frame of the container, during the active 
regime, we obtain the following equation of motion for 
the position y(t) of the intruder: 
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where u = 2-kJ , Ap = pi — p m , and u is a positive 
constant. 

We have assumed that the drag force is proportional 
to the density of the granular medium and to the velocity 
of the intruder. The constant vq may depend on the size 
and the shape of the intruder, and the effective viscosity 



of the medium while fluidized. In this model, the sand 
fluidizes when the term inside the parentheses, in equa- 
tion (1), is positive. This active regime occurs during 
the interval of time t when the acceleration, due to the 
moving noninertial reference frame, points upward and 
its magnitude is greater than gravity. As we mentioned 
before, this condition can be seen as an apparent grav- 
ity pointing upward, causing an increase in the relative 
interparticle motion. The time interval r, in one cycle 
of the oscillation, goes from the instant the granular sys- 
tem becomes active to the time when the passive regime 
starts, and is given by 
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At the beginning of this time interval r the velocity of 
the intruder relative to the noninertial reference frame is 
taken to be zero. 

Integrating twice equation (1) we calculate the dis- 
placement Ay T per cycle, during the time interval r. If 
we sum all the contributions per cycle in a time t much 
greater than the period T of the oscillations, the average 
vertical displacement y{t) of the intruder is obtained, 
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where a = lut/2, and we have neglected a transient term, 
assuming that vo(pi/p m ) ^> u>. The expression inside 
the square brackets is the average vertical velocity of the 
intruder. This velocity can be expressed as, 



where, 
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The model predicts the following characteristics of the 
migration of a large object in a granular medium when 
this fluidizes in part of the cycle: 

1. Heavy objects (pi > p m ) rise at constant velocity. 

2. Light objects (pi < p m ) sink at constant velocity. 

3. The dependence of the velocity on the ratio pij p m 
is linear and the slope is equal to the negative of 
the intercept, on the velocity axis. 

Near Ap = other contributions may become dom- 
inant, but it is assumed, for simplicity, that the model 
remains valid near p m ~ p. L . 

In figure 3 we can see that the experimental results 
agree with the predictions of the model. The slope is 
equal to the magnitude of the intercept, within two sig- 
nificant figures. If F(a, vq) = 12mm/ s, we get that the 
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value for the only free parameter is isq = 840s -1 . In 
the inset a plot of the velocity versus the radius of three 
large heavy spheres is shown. These measurements were 
made with three steel balls with different diameters, and 
the same density. This result suggests that the velocity 
v of the intruders depend linearly on the radius r of the 
spheres. This is consistent with a linear dependence of 
the drag force on the transversal area of the intruders, 
which implies that the retarding force is proportional to 
f 2 PmV- Further measurements arc needed to determine 
the dependence of the velocity on the size and shape of 
the intruder. 

No points are reported very close to Ap = 0. The 
reason is that, in such range, the velocity is very small 
and the motion becomes erratic. The same object can 
stay in one place for a long time and then drift slowly 
either upward or downward. 

Our results are significantly different to those reported 
by Mobius et al [13]. In their experiment all the intrud- 
ers rise, even if the density is less than of the granular 
medium. We believe that the reason for this is that they 
used tapping instead of sinusoidal oscillations to shake 
the container. We tested this hypothesis by shaking a 
container with sand, and using a light intruder. When 
the motion was approximately sinusoidal, the intruder 
would sink, and when we applied brief and distinct shocks 
the same intruder would rise to the top. A distinct shock, 
seen from the reference frame on the container, produces 
a large acceleration downward, and then a smaller accel- 
eration upward as the system relaxes back to its original 
position. This does not produce an effective reversal of 
the apparent gravity, during a sufficiently long time r. 
The mechanism involved in the migration of an intruder, 
may be different from the one found in a vibrated gran- 
ular matter, when tapping is used. 

We are presently studying in detail, the conditions for 
which reverse buoyancy can be observed. These condi- 
tions should give clues about the detailed mechanism of 
the cyclic fluidization and condensation of the granular 
medium subjected to vibrations. 

Although the assumptions arc crude, the predictions 
based on them work very well for our experiment. Wc 
see from the experimental results that density ratio can 
be a crucial parameter to describe the migration of large 
objects immersed in a medium of smaller particles. An 
important advantage of our model is that it gives a very 
simple explanation for the rising and sinking of large ob- 
jects in vibrated granular beds, in terms of the density 
ratio, for the case of reverse buoyancy. This experiment 
gives strong support to the idea that the oscillations can 
produce a series of active and passive regimes, where the 
granular medium can be considered a fluid, in part of the 
period of oscillation, and a solid in the remaining part. 

In conclusion we have provided a useful way of inter- 
preting the physical origin of reverse buoyancy through 
measurements of the velocity of a large object immersed 



in a vibrated granular bed, and a simple model that in- 
dicates the role of a buoyancy force and the drag force. 
With this model we are able to predict the velocity of 
the intruder as a function of its density. According to 
the proposed model the observed behavior is a conse- 
quence of the fluidization of the granular medium in only 
part of the cycle. A detailed understanding of how the 
granular medium is fluidizcd and the relative influence of 
various parameters is not yet available. It would be inter- 
esting to investigate the influence of ambient gas, and the 
effect of changing the physical properties of the grains. 
In particular, the presence of ambient gas may be neces- 
sary for the fluid-like behavior during the active phase. 
If this is so, then reverse buoyancy might be expected 
to be absent in a sufficiently hard vacuum. Indeed, in 
some situations, hard vacuum has been observed to re- 
sult in an analogous suppression of heaping of vertically 
vibrated granular media [18]. 
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FIG. 1. Sequence of pictures from a digital video of the 
vibrating granular system, with an intruder. The top three 
images are of an intruder of low density pi = 0.26gr /cm 3 , 
from lefth to right t = 0s, t = 3.83s and t = 5.27s. The 
bottom three pictures are of an intruder of high density 
Pi = 7.7gr/cm 3 , from left to right t = 0s, t = 0.77s and 
t — 1.10s. The rectangles indicate the position of the in- 
truder. 

FIG. 2. The measured vertical displacement for (a) heavy 
intruders and (b) light objects, as a function of time. The 
linear fit shows the agreement between the model and the 
experiment. 

FIG. 3. Vertical velocity of the intruder as a function of the 
density ratio, pij p m . The solid line represents a linear fit. The 
equation resulting from the fit shows the complete agreement 
beteen theory and experiment (see Eq.4) . The inset shows the 
dependence of the velocity of a heavy intruder on the radius 
of the sphere, for three steel balls of the same density. 
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This figure "figurel.jpg" is available in "jpg" format from: 
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